Rotation shields chaotic mixing regions from no-slip walls 
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We report on the decay of a passive scalar in chaotic mixing protocols where the wall of the vessel 
is rotated, or a net drift of fluid elements near the wall is induced at each period. As a result the 
fluid domain is divided into a central isolated chaotic region and a peripheral regular region. Scalar 
patterns obtained in experiments and simulations converge to a strange eigenmode and follow an 
exponential decay. This contrasts with previous experiments [Gouillart et a/., Phys. Rev. Lett. 
99, 114501 (2007)] with a chaotic region spanning the whole domain, where fixed walls constrained 
mixing to follow a slower algebraic decay. Using a linear analysis of the flow close to the wall, as well 
as numerical simulations of Lagrangian trajectories, we study the influence of the rotation velocity 
of the wall on the size of the chaotic region, the approach to its bounding separatrix, and the decay 
rate of the scalar. 



Chaotic advection is often presented [HE] as an effi- 
cient alternative to turbulence for the mixing of viscous 
fluids or delicate substances, as often encountered for ex- 
ample in the food and pharmaceutical industries. The re- 
peated stretching and folding of fluid particles, induced 
for instance by mechanical stirrers, transforms an het- 
erogeneity into a striated pattern with high gradients of 
the concentration field. This greatly enhances molecu- 
lar diffusion, which accelerates the homogenization pro- 
cess. In contrast to turbulence, the velocity field that 
causes chaotic advection is spatially smooth and coher- 
ent in time. Fluid particles may therefore reside for long 
times in regions where stretching is much lower than av- 
erage, before wandering to the rest of the domain. The 
global mixing rate is slowed down by such fluid particles 
whose motion is not ergodic on the typical timescales of 
a mixing experiment. For example, "sticky" KAM is- 
lands trap particles in their vicinity for long times [3]. 
Recent theory [4], numerical simulations [5], and experi- 
ments [6j [7J [8] have shown that the vicinity of a no-slip 
vessel wall is a region of poor stretching that can slow 
down the whole mixing process. 

In this Letter, we present a class of time-periodic 
protocols exhibiting chaotic advection, but where the 
chaotic region is shielded from the influence of the wall. 
We recover the exponential decay of the scalar variance 
(Fig.Qd)) and experimental concentration patterns that 
converge to a periodic pattern after a few periods (Fig.j2|, 
as predicted by theories based on statistics of stretch- 
ing [9j [TDJ [11] or eigenmode analysis of the advection- 
diffusion operator [I2j [13] . On the contrary, the conver- 
gence to an eigenmode pattern was never observed on 
experimental timescales for a chaotic region that extends 
to the walls [6] [7]. The key ingredient is that a regu- 
lar region separates the chaotic region from the wall, as 
in Fig. [TJb) and (c). Rotating the wall at a constant 
angular velocity is a simple way of creating such a regu- 
lar layer that insulates the chaotic region from the wall. 



In the present study, we investigate the case of a slowly 
moving wall and predict the size of the regular region, as 
well as the dynamics near the separatrix between chaotic 
and regular regions. We compare results of an analytical 
model with numerical simulations. Finally, we investi- 
gate the influence of the phase portrait on the dynamics 
of mixing and deduce an optimum for the rotation veloc- 
ity. 




FIG. 1: Chaotic mixing in a closed vessel, (a)-(b) Experi- 
ments: a rod moves periodically following (a) a figure-eight 
and (b) an epitrochoid path, (c) Numerical simulation of 
the figure-eight path with a rotating wall. Insets: Poincare 
section (stroboscopic map) obtained numerically for the cor- 
responding Stokes flow (dotted line: rod path; filled circles: 
fixed points of interest; solid line: associated separatrices) . 
(d) Evolution of the standard deviation a(C) of the concen- 
tration, rescaled by its value at Unit, the time at which the 
mixing pattern emerges, for (a)-(c). 
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FIG. 2: [Color online] Experimental evidence of an eigenmode 
for the epitrochoid protocol. Top: probability distribution 
function of the concentration. Bottom: rescaled concentra- 
tion field at successive periods, showing the converged peri- 
odic pattern. 

Let us first present the basic experimental and nu- 
merical observations. In all protocols, we consider two- 
dimensional flows with a single periodically-driven cylin- 
drical stirrer, and an outer cylindrical wall that is either 
fixed or rotated at a constant angular velocity. Experi- 
mentally, the rod gently stirs viscous sugar syrup follow- 
ing either a figure-eig ht path (Fig.[jja)) or an epitrochoid 
(Fig.[TJb)), that is a loop with a smaller inner loop. The 
fluid viscosity v = 5 x 10 -4 m 2 • s _1 together with rod 
diameter i — 16 mm and stirring velocity U = 2 cm • s _1 
yield a Reynolds number Re = Ui/v ~ 0.6, consistent 
with a Stokes-flow regime. A spot of low-diffusivity dye 
(India ink diluted in sugar syrup) is injected at the sur- 
face of the fluid, and we follow the evolution of the dye 
concentration field during the mixing process (see Fig.[T]). 
The concentration field is measured through the trans- 
parent bottom of the vessel using a 12-bit CCD camera 
at resolution 2000 x 2000. Since our experimental set-up 
does not allow rotation of the vessel wall, we resort to 
numerical simulations of Stokes flow [14] for integrating 
Lagrangian trajectories of the figure-eight protocol with a 
rotating wall. A scalar concentration field is obtained by 
coarse-graining the positions of a large number (4 x 10 6 ) 
of trajectories, all initialized inside the same small spot. 
Figure [TJc) shows a numerically-computed concentration 
field for the case of a rotating wall. 

In all cases, the scalar blob, initially released in the 
center of the vessel, is stretched and folded into thin fil- 
aments, resulting in a complicated pattern that expands 
with time and gradually fills the chaotic region. The ex- 



tent of the chaotic region can be seen in the numerically- 
computed Poincare sections in the insets of Fig. [lla-c). 
For the figure-eight protocol with fixed walls (FigTjlJa)), 
the chaotic region spans the entire vessel. By contrast, 
the other protocols force fluid particles to follow reg- 
ular closed trajectories in the vicinity of the wall, ei- 
ther because of the rotation of the wall itself (Fig. [TJc)), 
or because the epitrochoid protocol imposes a net ro- 
tational motion of the fluid elements close to the wall 
(Fig.[]Jb)). A regular (non-chaotic) unmixed region thus 
separates the walls from the chaotic region in Fig. [TJb— 
c). This is the key factor responsible for the different 
evolutions of the standard deviation cr(C) of the concen- 
tration field shown in Fig. []Jd). Statistics of the con- 
centration field are measured in a large rectangle inside 
the chaotic region. In the case of the figure-eight pro- 
tocol with fixed wall, the standard deviation decays al- 
gebraically, whereas it decays exponentially in the two 
other cases. For the epitrochoid protocol, we have exper- 
imental evidences that the concentration fields rapidly 
converges to an eigenmode (see Fig. [2|: the probability 
distribution of the concentration and the concentration 
pattern itself rapidly become invariant. This periodic 
pattern is associated with the slowest- decaying Floquet 
eigenmode of the ad vect ion- diffusion operator [12 . By 
contrast, we described in [6j [7] how no permanent pat- 
tern can be attained in the case of a fully chaotic region, 
because an unmixed pool of fluid close to the wall slowly 
leaks into the center of the fluid domain. 

In [6j [7] , we showed how the algebraic variance decay 
is related to the algebraic dynamics of the flow in the 
vicinity of a parabolic fixed point located at the wall. In 
the following, we will show how the new phase portrait 
obtained in the presence of a net rotation within the ves- 
sel is characterized by a regular region near the wall of 
typical width d. The convergence of fluid particles to this 
region is exponentially fast, and we recover an exponen- 
tial decrease of the variance. We focus on the figure-eight 
rod-stirring protocol with a moving wall, where we can 
tune the rotation rate of the wall independently from the 
rod motion. In the limit of slowly moving walls, incom- 
pressibility and the no-slip boundary condition suggest 
writing the flow near the wall as a simple map, 

= e + n + A(9)y, y = y- \A\e)y\ (1) 

to leading order in y. Here (6(t),y(t)) is the position 
of a fluid particle at the beginning of a time inter- 
val, (<9, y) = (6(t + T),y(t + T)) its new position after one 
stirring period T, and we have supposed that the fluid 
is incompressible. The angle is measured counterclock- 
wise along the wall, < y <C 1 measures the distance 
from the wall, and Q/T is the wall angular velocity. The 
circular container is assumed to have unit radius. We can 
regard A(0) as the contribution of the rod motion to the 
velocity field near the wall. As such, A(0) does not vary 
significantly when Q changes, as long as Q remains small 
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FIG. 3: (a) The distance d(t) from the edge of the central 
mixed region to the wall, for different values of CL. (b) The 
exponential decay rate /jl(Q) of d(t) and the distance 
separating the hyperbolic fixed point from the wall, mea- 
sured from simulations (solid lines) and from Eqs. (|3| and Q 
(dashed lines). 

(this has been confirmed numerically). A(9) is positive 
where the motion of the rod and the wall reinforce each 
other (i.e. in the left part of the vessel in Fig. [TJc)) and 
negative where there is opposition between the wall and 
the rod (right part of the vessel). In particular, there 
exists a fixed point of the map ([!]) (a period- 1 point of 
the flow) given by 

A'(0*) = O, y* = -CL/A(0*) (2) 

where the rotating wall balances exactly the flow induced 
by the moving rod. This fixed point is at a distance 

d*(n) = n/\A(e*)\ (3) 

from the wall and is pictured as a small filled disk in 
Fig.[TJc). The linearization matrix for the dynamics near 
the fixed point has eigenvalues 1 ± /i with 

fjL = Q ^-A"(0*)/2A(0*) (4) 

where the argument in the square root is non-negative 
since A(6*) < and A" (6*) > 0. This is a hyperbolic 
fixed point and the approach along its stable manifold 
(the separatrix plotted in Fig. [TJc)) is given by 

y(t)~I/(0) exp(-nt/T) (5) 

to first order in CL. The approach to the fixed point is ex- 
ponential, at a rate proportional to the speed of rotation 
of the wall. For comparison, it was shown in [6j [7] that 
the approach to the wall scales as \jt for CL = 0. Fig- 
ure[3^c) displays the evolution of d(i), the distance sep- 
arating the scalar blob from the fixed point, for simula- 
tions with different CL. One clearly observes the predicted 
exponential approach, with a slower rate for smaller CL 
(eventually leading to an algebraic decay for CL = 0). In 
Fig. J3^b) we show the results for d*(Ci) and /jl(CL) denned 
in Eqs. ([3| and Q, as measured in numerical simula- 
tions. We observe a linear evolution of d*(CL) and fJb(Cl) 
with CL that agrees quantitatively with Eqs. Q and Q. 



Even a slow rotation of the wall (CL = 0.3 corresponds 
to a whole rotation of the vessel in 21 periods of the 
rod motion) decreases the radius of the chaotic region by 
25%. Altogether, the agreement of the above analysis 
with numerical simulations is excellent. 

The different rates of approach are reflected in a differ- 
ent time-evolution for the scalar variance. To show this, 
we have simulated the mixing of a scalar blob for different 
wall velocities (CI varying from 0.05 to 0.4). The evolu- 
tion of the normalized standard deviation of the scalar 
field is measured inside a central rectangle, and plotted 
in Figs. Qb-d). Fig. |4^a) offers a simplified sketch of the 
observed evolution. For a few initial periods (t < tmit) 5 
the mixing pattern emerges and diffusion is not yet effi- 
cient. (Figs. |4^b-d) show a spurious initial 'spike' since 
the concentration is not measured in the entire domain.) 
For times t > t^a we observe a rapid exponential de- 
crease of cr(C) for all Ct. After a time t sep , we observe a 
transition to a slower exponential regime with a rate a 
that increases with Cl. The transition time t sep decreases 
when CL increases (see Figs.Qb-d)). Increasing CL further 
leads to the collapse of this second exponential regime 
on the continuation of the initial exponential regime (see 
Fig.gd)). 

We now propose a scenario that accounts for these 
observations. Between £ init and t sep , diffusion starts to 
smear out dye filaments that are elongated by the rod 
motion. The decay of the variance stems from the dis- 
tribution of stretching experienced by fluid particles in 
the bulk, as described in [9 j. Since the influence of the 
wall on the velocity field is weak in the region spanned 
by the rod's path, the same decay regime is observed for 
different values of CL. At the time t sep when some fila- 
ments arrive in the vicinity of the hyperbolic point and 
its separatrix (visible as a 'knee' in Fig. [3|, fluid par- 
ticles that have stayed most of the time in the core of 
the chaotic region have all been stretched significantly, 
and most of the initial variance corresponding to such 
fluid particles has been exhausted. However, there ex- 
ists outside the central mixing pattern a pool of poorly 
stretched fluid delineated by the separatrix. For t > t sep , 
the decay of the distance d(t) to this separatrix is gov- 
erned by the exponential dynamics of Eq. (J5|, and be- 
cause of area-preservation an unmixed strip of width 
A(t) ~ d(t) - d(t + T) ~ -Td(t) must enter the central 
mixing region along the unstable manifold of the fixed 
point (see the white 'tongue' in Fig. [TJc)). The strip is 
then thinned to the Batchelor length £ = \/tt/A, where k 
is the molecular diffusivity and A is an average stretch- 
ing rate inside the central mixing region [11] [15j [16] , and 
finally wiped out by diffusion. 

The fluctuations of the mixing pattern in Fig. [TJc) are 
dominated by the white strips injected most recently, as 
well as a few folds of strips injected previously — portions 
of the strips that are stretched but not folded reach £ 
faster. We therefore estimate that the total white area 
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FIG. 4: Evolution of the standard deviation <j(C). (a) Sketch 
of the decay of cr(C) according to the proposed scenario; (b)- 
(d) Actual decay measured in simulations of the figure-eight 
protocol, for three different wall velocities Q; the continuous 
(resp. dash) line shows an exponential decay with rate fi/2 
(resp. A/2). 

inside the bulk is the dominant contribution to the vari- 
ance, and that it is proportional to \d(t)\ ~ exp(— fit/T). 
We observe indeed a decay of the variance at a rate close 
to \i for small values of Vt (Fig. Qb-c)), consistent with 
a rate of mixing dominated by the rate of approach to 
the hyperbolic fixed point. The central mixing process is 
potentially more efficient (A > /i), but it is starved by the 
boundaries. When the rotation of the wall is increased 
so that [i becomes of the same order as A or greater, 
stretching in the bulk is as slow as that near the separa- 
trix (or slower), so that the decay rate saturates at the 
same value A as during the initial decay phase. 

In summary, we have proposed a design principle for 
mixing protocols that protect the chaotic mixing region 
from the slowdown due to no-slip walls. Creating a net 
rotation of the fluid in the vicinity of the wall results in 
a mixed phase space with a regular region that insulates 
the central chaotic region from the wall. In this case, 
we retrieve an exponential decay of scalar variance. For 
the case of a rotating vessel wall, rotating faster may 
increase the mixing efficiency. However, we have seen 
that the variance decay rate is always bounded by that 
of the bulk, and that when increasing Q the total mixing 
region is smaller, so the influence of the separatrix is felt 
earlier. An understanging of this mechanism is a first 
step towards the design of efficient mixers. Depending 
on the application, it may for instance be appropriate to 
rotate the wall in such a way that — A, to match 

the rate of stretching in the bulk to that near the outer 
hyperbolic point. 

From a practical point of view, we note that for the 
epitrochoid protocol the rod itself generates the net ro- 
tation of the fluid elements, which qualitatively induces 
similar effects as above, but without the need for a rotat- 



ing wall. Obviously, the biggest drawback of the design 
principle presented here — to insulate the chaotic region 
from the walls — is the presence of an unmixed region 
near the wall. But this may be unimportant in applica- 
tions where one can selectively sample the mixture, for 
instance in pharmaceuticals. In any case, there is no 
other known practical way to avoid a wall-induced slow- 
down of the variance decay. Finally, from a more fun- 
damental point of view it is intriguing that in all cases, 
including the epitrochoid, the decay rate inside the mix- 
ing region is always rather small, as compared for in- 
stance to that of a baker's map. The origin of this slow 
mixing — the distribution of hyperbolic points, parabolic 
points, etc. [jj |9j [TU1 E] — is a key issue to investigate in 
the future. 
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